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1.  Introduction 

! 

In  order  to  significantly  improve  the  ability  to  detect  underground  nuclear  explo¬ 
sions  using  seismic  measurements  and  to  minimize  the  biases  and  uncertainties  associ¬ 
ated  with  yield  estimations,  there  has  been  a  recent  trend  towards  the  use  of  extended- 
band  seismic  data  in  regional  distance  ranges.  This  data  holds  out  the  potential  of  higher 
resolving  power  than  the  traditional  teleseismic  data  which  has  been  used  over  the  years 
for  monitoring  underground  nuclear  testing  in  the  Soviet  Union,  however  the  regional 
discriminants  in  use  today  have  been  determined  largely  in  an  empirical  fashion  from 
data  collected  in  the  United  States  and  their  applicability  to  other  regions,  such  as  the 
Soviet  Union,  is  not  completely  obvious.  This  is  due  to  the  fact  that  we  do  not  under¬ 
stand,  from  a  theoretical  seismological  standpoint,  many  of  the  properties  of  seismic 
amvals  that  are  used  in  regional  discriminants,  P„.  S..,,  P„  and  L,.,  thus  making  it  difficult  to 
extrapolate  these  properties  to  a  different  setting. 

Confidence  in  our  discrimination  and  yield  estimation  methods  is,  to  a  large  extent, 
dependent  upon  our  understanding  of  P,„  S,„  P,.  and  L„  propagation  and  the  roles  played  by 
frequency  dependent  anelastic  attenuation  and  the  depth  dependence  of  structural  elastic 
parameters  in  the  crust,  at  the  Moho  and  in  the  upper  mantle.  Suitable  velocity  gradients 
in  the  vicinity  of  a  major  structural  discontinuity,  such  as  the  Moho.  can  have  large 
effects  on  the  associated  seismic  arrival,  P,„  and  low  velocity  zones  or  zones  of  "random¬ 
ized  laminations"  within  the  crust  can  act  as  veiy  efficient  waveguides  to  tiap  seismic 
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arrivals,  such  as  P^.  Simple  ray-based  modeling,  that  has  been  very  successful  at  telese- 
ismic  distances,  can  produce  mconect  and  misleading  results  when  applied  to  regional 
problems.  In  addition,  we  have  found  that  even  "complete"  synthesis  methods  can 
mislead  us  if  we  fail  to  represent  the  fundamental  characteristics  of  hypothetical  struc¬ 
tural  models  correctly. 

In  this  paper  we  report  on  the  initial  results  of  a  long  term  study  which  is  aimed  at 
obtaining  understanding  of  the  fundamental  proccs.ses  involved  m  legional  seismic  wave 
propagation  and  what  the  data  that  we  observe  tell  us  about  the  nature  of  the  eailh’s  ci  ust 
and  upper  mantle.  With  this  knowledge  we  w'lll  be  able  to  more  confidently  resolve 
detailed  source  characteristics  and  make  it  possible  to  significantly  improve  detection  and 
yield  estimation  capabilities. 

Our  basic  objective  is  to  develop  an  inversion  algorithm  which  will  directly  com¬ 
pare  broadband  regional  data  with  complete  synthetic  seismograms  to  infer  detailed 
structure  and  source  properties.  This  objective  is  ambitious  and  touches  on  most  areas  of 
seismology,  both  observational  and  theoretical.  As  a  stalling  point,  we  began  by  using 
the  results  of  other  researchers,  such  as  Gomberg  and  Masters,'  who  have  successfully 
inverted  for  crust  and  upper  mantle  structural  parameters  by  directly  comparing  complete 
kK’ked  mode  synthetic  seismograms  with  the  observed  data  in  the  lime  domain.  The  kov 
to  this  inversion  is  the  use  of  synthetic  differential  seismograms  which  describe  the 
linearized  relationship  between  model  parameters  and  the  resulting  synthetic  seismo¬ 
grams.  These  provide  the  T’lechet  derivatives  that  are  necessary  m  the  inversion  and  it  is 
important  that  they  be  computed  accurately  and  efticiently. 

Prev  ious  inversion  efforts  of  this  type  have  all  been  limited  to  low  frequencies  (less 
than  0.2  Hz)  and  have  directly  compared  synthetics  with  data  m  the  time  domain.  Our 

'  Gomberg,  J  and  Masteis,  T,  088.  Wavefoiiii  nuxtelling  using  Ux;ked-nh\le  svmheiic  and 
ditfeieniuii  seismograms  application  to  deicrmiiialion  of  ilie  stiuctuic  of  Mcxuo,  Geoplivsical 
lou.'Tial  of  the  Royal  Astronomical  Society.  V  94,  p  |•')'^-2I.8 
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objective  is  to  push  this  inversion  to  the  highest  frequencies  that  will  produce  useful 
results  and  wc  would  certainly  hope  to  get  to  at  least  1  hz.  In  order  to  accomplish  this  it  is 
necessary  to  solve  the  following  key  problems. 

•  The  development  of  an  accurate  and  efiicicnt  algorithm  for  computing  differential 
seismograms  which  are  necessary  for  the  inversion  procedure. 

•  The  development  of  better  methods  for  incorporating  anelastic  attenuation  into 
modal  synthesis  computations. 

•  The  determination  of  starting  solutions  that  will  produce  synthetic  seismograms 
which  have  the  same  general  characteristics  that  we  see  in  the  data. 

In  section  2  we  discuss  our  results  in  developing  an  exact  differential  seismogram 
program.  This  was  necessary  because  the  traditional  first  order  penurbation  method  for 
computing  differential  seismograms  proved  to  be  inadequate  for  our  problem.  In  section 
3  we  talk  about  how  we  have  dealt  with  the  problems  related  with  anelastic  attenuation 
and  in  section  4  we  show  results  that  give  us  a  powerful  new  way  of  looking  at  the 
earth's  crust  which  provides  us  with  an  easy  method  for  modeling  many  of  the  charac¬ 
teristics  that  we  see  m  broad  band  regional  data. 

2.  Computation  of  Differential  Seismograms  at  High  Frequencies 

The  standard  method  for  computing  differential  seismograms  makes  use  of  first 
order  peilurbation  (F-OP)  theory  and  is  described  by  Takeuchi  and  Saito.  In  this  tech¬ 
nique  Rayleigh’s  principle  is  used  to  find  expressions  for  the  effects  of  small  changes  of 
the  elastic  properties  on  the  phase  velocities  of  normal  modes.  These  expressions  relate 
the  partial  derivatives  of  the  eigenvalues  to  depth  integrals  of  the  un-perturbed  eigen¬ 
functions  and  make  it  possible  to  compute  eigenvalue  derivatives  implicitly  without 
direct  application  of  the  chain  rule  tliiough  the  entire  sequence  of  computations.  The 

-  Takeuclii,  H  and  Saito,  M  ,  1972.  Seismic  Suiface  Waves  Methods  in  Computational  Phy¬ 
sics,  V,  1  I ,  p  21  '-294,  ed  Holt,  B  .  Acadenne  I'less.  New  York 
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resulting  formalism  makes  it  possible  to  compute  eigenvalue  derivatives  efficiently  and 
accurately  and  this  method  has  been  used  extensively  throughout  the  seismological  com¬ 
munity.  However,  the  variational  principle  does  not  produce  eigenfunction  derivatives  so 
use  of  this  method  implies  that  the  eigenfunction  derivatives  will  be  neglected  in  the  final 
differential  seismograms. 

In  cylindrical  coordinates  we  can  express  the  frequency  dependent  P-SV  displace¬ 
ment  vector,  Ur„  for  a  single  normal  mcxle,  n,  as  follows. 


u,(r,.0,,z,)  =  k„£  [L(k„.m)]  |£(ri,z,)l  E ,(n.z,)P(k„.in,i,,0,)  +  £:(ii,z,)B(k„.ni,i,,tV);  i 

rn  V 


(2.1) 


where. 


ni  IS  the  azimuthal  order  number, 
r,  IS  the  radial  distance  to  the  receiver, 
e,  IS  the  azimuth  to  the  receiver, 
z,  IS  the  depth  of  the  receiver, 
z,  IS  the  depth  of  the  source, 

k  IS  the  frequency  dependent  eigenwavenumber  for  mode  n. 

[I]  IS  the  frequency  dependent  four  competent  source  jump  vector, 

[£  ]  IS  the  frequency  and  depth  dependent  four  component  eigenfunction  vector, 
P,B  are  the  P  and  B  vector  cylindrical  harmonic  components 


By  simple  application  of  the  chain  rule  we  can  express  a  differential  displacement, 
c)u  y()v,  where  v  is  a  model  parameter,  as  follows. 


du.yuv  =(fJk./avj  ju:,/k:,  +  k,.V  |(dlL]/dkj|£  I  i/iiP  +  Z:  Bj  +(L11£  I  i£:(dPAJk  )  +  £  (dBMkj! 
-t-k.V  [Lii'a|£  l/dv)  'A.  ,P  +  A'_B  +  \L\\E  \  Udli./dw)P  +  QH  Eh  )B:| 


When  using  FOP  theory  to  compute  differential  seismograms,  only  the  first  term  in  equa¬ 
tion  (2.2)  which  depends  on  the  eigenvalue  derivative,  (.dkjfhn,  is  used  and  the  remainiag 
terms  in  (2.2)  which  depend  on  the  eigenfunction  derivatives.  !d|/  i/dvi.idHi/dv)  and 
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are  igi.ored. 

Our  initial  implementation  of  a  differential  seismogram  synthesis  program  used  the 
standard  FOP  approach.  Since  vve  were  planning  on  using  the  differential  seismograms  at 
relatively  high  frequencies,  we  fell  that  it  would  he  prudent  to  carefully  check  the  vali¬ 
dity  of  the  approximation  that  we  were  using.  We  did  this  by  computing  a  set  of  difter- 
ence  seismograms  where  we  took  an  initial  structural  model  and  computed  locked  mode 
seismograms,  made  a  small  change  in  a  single  model  parameter  and  computed  new 
locked  mode  synthetic  seismograms,  differenced  the  new  seismograms  with  the  original 
seismograms  and  repeated  this  process  for  the  other  model  parameters.  In  this  way  we 
were  able  to  obtain  approximations  to  the  exact  first  order  Taylor  senes  terms  which 
included  all  terms  in  equation  (2.2)  and  which  were  numerically  accurate  as  long  as  the 
parameter  changes  were  small  enough  so  that  the  first  order  term  in  the  Taylor  series 
dominated  the  higher  order  terms. 

Figure  I  shows  the  results  of  a  comparison  of  analytic  differential  seismograms, 
labeled  H4P().s,  which  neglected  the  eigenfunction  derivatives  .  versus  difference  seismo¬ 
grams.  labeled  SdPO.i.  These  seismograms  weie  computed  for  a  simple  layer  over  a  half 
space  structure  representing  a  crustal  layer  over  the  upper  mantle  and  the  difference  and 
differential  seismograms  are  with  respect  to  the  f’-wave  velocity  in  the  crust  layer.  TTe 
bandwidth  of  these  seismograms  is  0-2  Iv.  and  the  source  receiver  distance  is  500  km. 
From  this  tigiire  it  is  obvious  that  the  analytical  differential  seismogram  using  FOP 
theory  is  m  error.  The  time  window  shown  in  figure  I  repiesents  the  early  pan  of  the  P, 
wave  tram  and  when  we  looked  at  ,S-wave  differentials  later  during  the  Rayleigh  wave 
w'e  got  good  agreement.  'Hus  shows  that  F'OP  theory  is  not  adequate  for  computing  dif¬ 
ferential  seismograms  at  frequencies  aiound  I  hz. 

In  ordei  to  convince  ourselves  that  the  problem  was  associated  with  neglecting  the 
eigenfunction  derivatives,  we  conducted  an  experiment  wheie  we  replaced  the  eigenfunc¬ 
tions  for  the  perturbed  seismograms  befoie  differencing  with  those  of  the  unperturbed 


Time  (sec) 


Figure  1.  Difference  (S4P05)  vs.  differential  (H4P05)  seismograms; 
eigenvalues  only. 
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Figure  2.  DiffcieiKC  (S4P().‘>)  vs  dilfeiential  (S.-XKi.s)  seisnutgiams: 
eigenvalues  and  eigenfunctions. 


seismograms  and  recomputed  ditTerence  seismograms  where  the  eigenfunctions  were 
forced  to  remain  constant.  In  this  case  comparisons  of  our  modiHed  difference  seismo¬ 
grams  and  the  analytical  FOP  differential  seismograms  were  very  good,  indicating  that 
the  problem  was  associated  with  changes  in  the  eigenfunctions. 

We  then  developed  a  computer  code  which  computes  differentiai  seismograms 
analytically  without  neglecting  the  eigenfunction  derivatives.  This  proved  to  be  a  tedious 
and  difficult  task,  however  we  were  able  to  develop  a  program  w'hich  is  accurate  and 
relatively  efficient,  ceitamly  when  compared  to  the  differencing  approach.  Our  method  in 
developing  this  program  was  straightforward; 

•  We  abandoned  using  a  variational  principle,  which  involves  depth  integrals  of  the 
eigenfunctions,  for  computing  the  eigenvalue  derivatives  and  replaced  this  approach 
with  explicit  derivative  computations  using  the  chain  rule.  This  was  a  tedious  and 
laborious  process  which  involved  carrying  derivative  computations  through  the 
entire  chain  of  algebraic  operations,  however  we  found  that,  after  some  reairange- 
ment  and  algebraic  simplification,  the  resulting  numerical  algorithm  was  relatively 
efficient. 

•  We  then  checked  that  the  eigenvalue  derivatives  matched  those  from  the  variational 
computations  and  we  also  checked  against  numerical  eigenvalue  denvatives  from 
our  difference  seismograms. 

•  .Since  the  derivative  computations  had  been  cairied  through  all  of  the  intermediate 
steps,  it  was  then  relatively  straightforward  to  extend  these  computations  to  produce 
eigenfunction  derivatives. 

•  I  he  modal  summation  expressions  were  then  modified  to  conform  w'lth  equation 
(2.2)  to  produce  exact  differential  seismograms. 

Figure  2  shows  a  comparison  of  our  exact  differential  seismogram' ,  labeled  SAIO.*). 
with  the  difference  seismograms,  labeled  .S2P0.‘i.  it  is  obvious  from  this  figure  that  the 
comparison  is  very  good  and  it  is  likely  that  the  residual  error  is  due  to  the  numerical 


approximation  implicit  in  the  difference  seismogiam.  An  example  showing  exact  dif¬ 
ferential  seismograms  for  a  realistic  regional  situation  with  a  complex  structural  model 
can  be  seen  in  figures  3  -  6.  The  stiuctural  P  and  S  velocities  are  shown  in  figure  3  for  a 
70  layer  model.  As  can  be  seen,  the  velocity  vs.  depth  profiles  have  a  random  component 
which  we  will  discuss  in  detail  in  section  4.  This  structure  represents  a  starting  estimate 
for  the  region  around  the  Soviet  nuclear  test  site  at  Semipalatinsk,  KSSR.  Figure  4  shows 
differential  seismograms  with  respect  to  P-wave  velocities  for  each  layer  in  the  model. 
The  70  differential  seismograms  are  plotted  one  above  the  other  as  a  function  of  layer 
index  and  the  original  synthetic  seismogram  is  also  shown  at  the  top  of  the  figure  for 
reference.  Figure  5  is  the  same  as  figure  4  except  that  the  S-wave  differential  seismo¬ 
grams  are  shown.  Figure  6  is  a  repeat  of  figure  4  except  the  vertical  axis  is  layer  index 
instead  of  depth  (for  comparison  with  figures  3  and  4. 

Figures  3  and  4  give  us  much  information  about  the  nature  of  the  regional  seismic 
wave  propagation  for  this  example.  The  source-receiver  distance  is  254  km  and  we  can 
see  how  different  regions  within  the  crust  effect  the  resulting  seismograms.  The  L,  coda 
IS  primarily  controlled  by  the  very  near  surface  part  of  the  ciust  although  the  initial  por¬ 
tion  of  the  aiTival  responds  to  the  entire  crust.  The  P,.  coda  is  also  strongly  effected  by 
near  surface  stiucture  and  we  can  clearly  see  the  strong  direct  Moho  rellection.  P„P.  as  a 
vertical  streak  down  to  layer  number  35.  Both  P„  and  S„  can  also  be  seen  in  the  differen¬ 
tial  seismograms  at  the  bottom  of  the  model. 

3.  Incorporation  of  Anelastic  Attenuation  in  Locked  Mode  Synthetic  Seismograms 

The  standard  method  for  incorporating  anelastic  attenuation  m  modal  summation 
seismogram  synthesis  methods  has  always  involved  the  use  of  FOP  theory.  .As  we  have 
seen  with  the  differential  .eismograms,  there  is  reason  to  question  ihe  aecuraev  of  F\')P 
theory  in  accounting  for  Q  effects.  Other  reseaichers  have  encounteied  ihis  problem^ 

’’  Day,  S  ,  McLaughlin,  K  ,  Shkoller,  B  and  Stevens,  J  ,  1989,  Poieniial  eirois  m  locked  mode 
synthetics  for  anelastic  earth  models,  Geopiiysical  Reseaich  letters,  v  16,  p,  2rrt-20o 
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i'lgure  6.  Stiuctuic  from  fig.  3  vs  layor  number 
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and,  as  with  the  differential  seismograms,  the  problem  seems  to  become  more  pro¬ 
nounced  at  higher  frequencies  and  higher  phase  velocities. 

In  our  early  efforts  to  address  this  problem,  before  we  had  solved  the  differential 
seismogram  problem,  we  developed  a  method  for  locating  the  complex  eigenvalues 
exactly  along  w  ith  a  complex  version  of  the  eigenfunction  and  modal  summation  codes. 
These  programs  work  very  well,  even  in  situations  where  the  Q  is  low,  although  the  ver¬ 
sions  of  the  programs  that  we  initially  developed  are  not  very  efficient.  Our  success  with 
the  differential  seismograms  provides  another  track  for  solving  the  modal  Q  problem. 

Q  corrections  using  FOP  theory  involve  the  compulation  of  eigenvalue  shifts  result¬ 
ing  from  complex  elastic  parameter  shifts  due  to  the  intrinsic  attenuation.  This  normally 
results  in  purely  imaginaiy  shifts  of  the  eigenwavenumbers'^  which  are  then  represented 
in  the  modal  summation  as  frequency  dependent  decaying  exponential  temis.  Thus  nor¬ 
mal  modal  Q  corrections  ignore  the  eigenfunction  shifts  m  the  same  manner  as  FOP  dif¬ 
ferential  seismograms.  In  high  frequency  -  high  phase  velocity  situations,  the  resulting  Q 
collections  are  in  error  even  for  very  high  Q  values. 

We  have  started  the  development  of  an  "exact  first  order"  modal  Q  collection  algo¬ 
rithm  based  upon  the  eigenfunction  derivative  capabilities  that  we  developed  for  the  dif¬ 
ferential  seismograms.  Tins  will  provide  an  efficient  means  for  computing  modal  Q 
corrections  that  will  always  be  accurate,  as  long  as  the  Q  values  are  sufficiently  high  so 
that  the  first  order  expansion  is  valid.  Initial  results  from  this  work  indicate  that  we  are 
close  to  achieving  this  goal. 

4.  Velocity  Randomization  and  Its  Effects  on  Regional  Synthetic  Seismograms 

When  investigating  an  extended  band  regional  seismogram,  one  is  souck  by  the 
large  amount  of  information  that  seems  to  be  contained  w  ithin  the  seismic  signal.  One 

^  Tlie  sliitts  are  purely  miapmaty  as  luiic  ,is  IJ  lelated  dispeisioii  effects  in  ihe  real  parts  of  the 
elastic  wave  vehxtities  are  iieplecied 
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does  not  see  simple  ray  arrivals  with  well  defined  coherent  wavelets,  but  instead,  an  ivals 
characterized  by  complex  wavetrains  with  onset  limes  that  are  often  difficult  to  deter¬ 
mine.  Typical  regional  .seismograms  that  were  recorded  near  the  Soviet  nuclear  test  site 
at  Semipalatinsk  are  shown  in  figures  7  and  8.  I'lgure  7  shows  vertical  components  at 
four  stations  from  the  Soviet  JVE  nuclear  test  as  a  function  of  receiver  range  from  the 
shot  site.  These  are  all  broad  band  seismograms  (1  -  .SO  hz).  The  stations  at  170  km  and 
2.53  km  were  the  temporarily  re-occupied  NRDC  sites  of  Karasu  (KSl  )  and  Karkaralmsk 
(KKL)  and  the  stations  at  1350  km  and  1529  km  were  the  IRIS  stations  at  Chusal  (CHS) 
and  Arti  (ARU).  Figure  8  was  taken  from  an  American  Cleophysical  Ihiion  poster  ses¬ 
sion^  and  shows  the  CHS  and  ARLI  data  after  applying  a  0.8  -  2  hz  passband  filter,  as 
well  as  additional  data  from  a  Peaceful  Nuclear  Explosion  (PNE)  which  took  place  in 
northwestern  USSR. 

These  data  are  typical  of  regional  explosion  events  m  the  shield  region  of  central 
and  northern  asta  and  from  these  data  we  can  make  the  following  observations. 

•  The  P„  aiTival  is  energetic  and  well  developed  at  all  distances  above  250  km  whore 
it  first  emerges  from  the  P,  wavetrain. 

•  L,.  IS  clearly  defined  for  distances  less  than  2(XX)  km  but  disappears  above  25lX)  km. 
The  L,  to  P„  amplitudes  stay  about  one. 

•  Apparent  Q  values  are  relatively  high.  High  frequency  energy  propagates 
efficiently. 

Most  researchers  attribute  the  incoherence  (9' regional  wavetrains  to  "lateral  scatter¬ 
ing"  which  IS  normally  intended  to  mean  some  three  dimensional  distribution  of  struc¬ 
tural  inhomogeneities.  If  we  assume  a  simple  uniform  distribution  of  seattercrs  m  space, 
at  a  variety  of  wavelengths,  then  the  scattering  of  a  coherent  wavefront  propagating 

ruveii,  H  ,  Itedlin,  M  ,  Rvigvi.  I  .  X'cinoii.  I- .  Kappiis,  M  .r'liaxt'/.,  I)  ,iikt  .VsIli.  P  ,  Ktiium 
al  stisaiic  observations  of  luicle.ii  e.xplosrons  inside  ilie  .So\  lel  I  ’nion,  I'ostei  I'lesenUkl  ai  ihc 
Aineiiean  Geophysical  Union  ineeiinp.  Deceiiilsei  1988.  abstiacl  appeared  in  KG.S.  I'l.ins  .Ainei 
Geophys  Union,  v  69,  p  1321 
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through  the  medium  will  result  m  an  attenuation  effect  of  the  coherent  wavefront  along 
with  incoherent  scattered  energy  that  follows  the  coherent  wavefront  in  a  coda.  The  net 
effect  is  that  the  energy  gets  "smeared  out"  spatially  so  that  arrivals  are  weaker  than  they 
would  have  been  m  a  smooth  medium. 

A  weak  unifomi  scattering  medium  will  result  in  high  apparent  Q  values  with  short 
and  weak  attendant  coda  and  a  strong  uniform  scattering  medium  will  result  m  low 
apparent  Q  values  with  long  and  energetic  coda.  The  data  we  obseiwe  in  the  Soviet  Union 
indicates  relatively  high  Q  values,  which  indicates  a  weak  scattering  medium,  yet  the 
coda  are  strong,  which  indicates  a  strongly  scattenng  medium.  A  way  of  resolving  this 
dilemma  is  to  look  for  other  scattering  mechanisms  which  allow  for  efficient  propagation 
of  seismic  energy  while  scattering  the  coherent  arrivals  strongly  to  produce  the  coda  we 
obseiwe. 

We  have  been  investigating  such  a  scattering  mechanism  which  assumes  that  the 
structural  inhomogeneities  are  anisoiroptc,  i.e.  the  scale  length  of  the  inhomogeneities  is 
different  m  the  veHical  direction  than  it  is  in  the  horizontal  direction.  Tlie  simplest  such 
anisotropic  scatterers  to  model  are  those  which  are  uniform  m  the  horizontal  direction 
and  arbitrarily  inhomogeneous  in  the  vertical  direction.  Of  course,  this  modeling  capabil¬ 
ity  has  been  in  existence  for  some  time  now,  however  researchers  using  laterally  homo¬ 
geneous  synthetic  seismograms  have  always  used  smooth  or  what  we  would  call  "large 
scale  blocky"  structural  models  which  will  not  produce  the  types  of  scattering  that  we 
would  expect  m  a  "random"  medium. 

The  results  from  a  traditional  smooth  model  are  shown  m  figure  10,  which  displays 
vertical  component  synthetic  seismogiams  as  a  function  of  source-receiver  distance  for 
the  structural  model  given  in  figures  9a  and  9b.  The  amplitude  scales  arc  adjusted  by  a  r- 
factor.  'I'he  structural  model  represents  the  eastern  Kazakhstan  region  and  we  represent 
velocity  gradients  with  homogeneous  lavei  approximations.  The  velocity  gradients 
helped  considerably  in  boosting  the  iclative  amplitude  of  P  over  what  it  was  without  the 
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Figure  10.  Syntlictic  vertical  components  for  the  base  Kazakh  model  in  fig. 
da-9b. 
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Syntlielic  vertical  components  toi  the  randomized  model  in  I'lg. 
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gradients,  however  the  seismograms  are  simple  when  compared  with  the  data. 

An  example  of  synthetic  seismograms  produced  fiom  a  vertically  "randomized" 
structure  is  shown  m  figure  1 1.  I'he  amplitude  scales  are  the  same  as  those  for  hguie  10. 
The  structure  model  is  shown  m  figures  rte-‘)d  and  was  determined  by  using  the  smooth 
model  in  figures  ^a-Ob  and  applying  a  landom  tluctuation  with  depth  dependent  layer 
thickness  and  variance.  "We  decided  to  use  100  meter  layer  thicknesses  for  the  upper  2 
km  and  1000  meter  layer  thicknesses  elsewhere.  We  also  used  high  variances  near  the 
stiiface  and  decreased  the  variance  with  depth.  Our  intent  was  to  crudely  match  the  sort 
of  near  suiface  velcx;ity  variances  seen  m  well  logging  measurements. 

A  comparison  of  figure  1 1  with  figure  10  shows  many  interesting  characteristics  of 
anisotropic  scattering  which  are  given  below. 

•  Strong  and  persistent  ccxla  are  generated  throughout  the  seismograms  for  the  verti¬ 
cally  randomized  structure.  The  later  arriving  phases.  P.,.  S,.  and  L,.  although 
coherent  and  easy  to  identify  in  the  smooth  stmeture  are  incoherent  with  ambiguous 
onset  times  in  the  vertically  randomized  structure. 

•  The  seismograms  for  the  vertically  randomized  structure  are  much  more  energetic 
throughout  the  duration  than  those  of  the  smooth  structure.  Not  only  has  strong 
ccxia  been  generated  in  the  "dead"  regions  of  the  seismogram,  but  the  amplitudes  of 
all  of  the  initial  airivals  have  been  preseiwed  as  well. 

The.se  results  are  significant  m  that  they  show  how  anisotropic  scattering  is  funda¬ 
mentally  different  from  isotropic  scatteiiiig  Instead  of  attenuating  and  spreading  out  the 
signals  as  a  unifomi  scattering  material  would  do,  a  vertically  randomized  medium 
focuses  the  seismic  energy  while  scatteiing  it  at  the  same  time.  The  horizontal  lamina¬ 
tions  introduce  numerous  low  velocity  zones  with  many  horizontal  layer  interfaces  that 
tend  to  flap  the  seismic  energy  similar  to  the  way  a  fiber  ciptic  cable  traps  light  This 
energy  is  free  to  propagate  horizontally,  however  it  is  inhibited  from  propagating  verti¬ 
cally  and  thus  energy  that  would  nomially  pre^pagate  away  through  the  bottom  of  the 
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structure  is  kept  concentrated  in  the  crust  and  upper  mantle.  dTie  introduction  ot'  lionzon- 
tal  laminations  within  a  region  effectively  inlioduce  a  negative  Q  eftect  since  they  will 
tiap  energy  within  the  laminations  and  thus  oveicome  the  normal  three  dimensional 
geometnc  spreading  effect.  The  focusing  associated  with  sliuciuial  laminaiions  is  dif¬ 
ferent  from  focusing  that  is  normally  associated  with  other  three  dimensional  lenses 
■Most  structural  lenses  focus  seismic  energy  into  small  and  well  delincd  regions,  however 
structural  laminations  channel  seismic  energy  into  broad  horizontal  sheets  so  that  the 
focusing  effect  can  be  seen  over  large  distance  ranges. 

The  ip.'roduction  of  veitical  randomization  into  smooth  stiuclural  models  uill  play 
an  important  role  in  explaining  regional  seismograms.  This  gives  us  an  c.xplainable. 
plausible  and  implementable  method  for  modeling  the  features  that  we  see  in  the  data, 
such  as  the  apparent  contradiction  of  high  Q  values  and  strongly  scattered  coda.  We  think 
that  this  will  also  stiongly  effect  how  Q  estimates  are  made  and.  ultimately,  yield  estima¬ 
tions.  We  can  see  from  these  examples  that  vertical  randomization  can  have  very  pro¬ 
nounced  effects  on  seismic  energy  levels  that  could  otherwise  be  interpreted  incoriectly. 
We  were  able  to  pump  up  the  P,.  coda  arbitrarily  by  putting  large  landom  tUictualions  at 
the  depth  where,  from  the  differential  seismograms,  we  knew  the  P.  energy  navels  We 
are  currently  testing  new  models  m  which  we  have  put  regions  of  high  lamination  at  the 
Moho  depth  in  order  to  pump  up  the  P„  portion  of  the  seismogram. 

In  figure  12  we  show  a  comparison  of  real  data  with  a  synthetic  seismogiam  that  we 
computed  using  the  randomized  structure  given  in  figures  9c-^)d.  This  shows  the  \ertical 
component  at  KKL  which  was  about  250  km  fimm  the  shot  and  both  the  data  and  the  syn¬ 
thetic  have  been  filtered  so  that  the  passband  of  comparison  is  about  1  to  2  hz  I  his  is  one 
of  the  best  real-synthetic  matches  of  regional  data  in  this  fiequencv  lange  that  we  have 
seen  and  we  think  that  this  type  of  qualitative  match  is  necessaiy  hefoie  a  formal  iinei- 
sion  procedure  can  be  reasonably  expected  to  produce  meaningful  results 


l-'iguic  12.  Comparison  of  real  data  (KKUYl-')  and  synthetic  (SYNOl)  for  the  JVli  at  station  KKl 
using  the  stiuctuie  in  tig.  9c-9d. 


5.  Conclusions 


In  order  to  improve  our  abilities  to  discriminate  low-yield  nuci  ir  explosions  and  to 
obtain  accurate  yield  estimates,  we  have  undertaken  a  study  to  inter  detailed  source  and 
structure  parameters  by  direct  inversion  of  broad  band  regional  seismic  data  using 
laterally  homogeneous  forward  modeling  methods.  We  have  identified  three  key  prob¬ 
lems  which  must  be  solved  in  order  to  accomplish  our  objective:  the  development  of  an 
accurate  and  efficient  algorithm  for  computing  differential  seismograms  which  are  neces- 
sai7  for  the  inversion  procedure,  the  development  of  belter  methods  lor  incorporating 
anelastic  attenuation  into  modal  synthesis  computations,  and  the  determinatk'n  of  starting 
solutions  that  will  produce  synthetic  seismograms  which  have  the  same  general  charac- 
tenstics  that  we  see  in  the  data. 

We  have  made  significant  progiess  in  each  of  these  problems.  We  determined  that 
the  normal  first  order  perturbation  theoiy  method  for  computing  differential  seismograms 
was  inadequate  for  higher  frequencies  and  phase  .  .nes  and  we  have  developed  and 
tested  an  exact  analytic  differentia'  .ek^mogiam  program.  We  also  have  developed  an 
e.xact  complex  pole  based  m  'thod  for  computing  modal  seismogiams  which  avoids  ihe 
use  of  first  older  perturbation  theory  for  incoipoiatmg  anelastic  attenuation  We  are  in  ihc 
process  of  developing  an  exact  first  order  Q  collection  ulong  the  lines  ot  the  exact  dif¬ 
ferential  seismogram  program  whicl  we  hope  will  piovide  an  efticient  means  for  incor¬ 
porating  Q  effects. 

Our  most  significant  discoveiy  is  that  anisotropic  scattering  pioduces  the  types  of 
effects  that  can  be  seen  in  the  observed  data  We  consider  this  to  be  oui  mosi  important 
result  to  date  because  it  f  as  proven  to  be  difticuli  te.  produce  sxnthelic  seismogiams  that 
match  the  gro.-.s  characterist  cs  of  broad  band  legional  data  Main  ivseaichers  attiibule 
this  difficulty  to  complex  and  mdelermmaie  lateral  scattering  processes  if  ihex  aic  right, 
then  It  IS  unlikely  that  we  will  iindeistand  ihe  nature  of  regional  nave  propagation  and 
we  will  be  forced  to  resort  almost  entirely  to  empirical  methods  for  doing  m-counti\ 
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yield  estimation.  We  have  discovered  that  an  explainable,  plausible  and  implementable 
method  for  modeling  the  features  that  we  see  in  the  data  is  to  use  vertically  randomized 
structural  models  in  laterally  homogeneous  modeling  codes. 

We  think  that  the  idea  of  horizontal  laminations  within  the  earth  is  very  plausible 
and  we  have  demonstrated  that  the  consequences  are  sigmficant.  We  know  from  well  log 
measurements  that  the  near  surface  earth  structure  looks  like  a  stochastic  process  with 
depth.  We  can  only  speculate  about  the  detailed  nature  of  the  structure  at  depth,  however 
it  is  unreasonable  to  rule  out  the  possibility  of  laminated  structures  deeper  in  the  earth.  A 
reasonable  hypothesis  for  oansition  zones,  like  the  Moho,  is  a  region  where  the  matenal 
on  one  side  gradually  "feathers"  into  the  material  on  t^he  other  side  to  produce  a  lamina¬ 
tion  zone  with  many  thin  layers  of  different  matenal  properties.  Such  a  zone  could  look 
like  a  velocity  gradient  in  travel  time  studies  and  for  longer  wave  length  waves,  however 
at  higher  frequencies  the  seismic  energy  would  be  efficiently  trapped  within  the  lamina¬ 
tions. 

With  the  results  from  our  studies  to  date  we  feel  that  we  are  in  a  good  position  to 
continue  our  efforts  toward  the  development  of  full  wave  regional  inversion  capabilities. 
We  aie  continuing  our  work  in  this  area  and  we  feel  confident  that  we  will  realize  our 
objectives. 
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